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The massless vector Schwinger model with Nf = 1, 2, 3, 4 number of flavors is studied on the lattice using the 
overlap formalism. A full Monte Carlo simulation yields values for the bilinear fermion condensate that are in 
agreement with the exact solution of the continuum Schwinger model with an added Thirring interaction. 
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The solution to the massless vector Schwinger 
model is known in the continuum. There the 
role of topology is clear: In the single flavor case, 
(\J/\J/) = if the gauge field has topological charge 
? ^ ±1 and ^ if q = ±1. No previous 

lattice regularization of chiral fermions could re- 
produce this result but the overlap formalism [1] 
can. Here we present the results of a MC simu- 
lation of the massless vector Schwinger model on 
the lattice for Nf = 1, 2, 3, 4 flavors in the overlap 
formalism. Preliminary results were reported in 
[2]- 

The continuum action in Euclidean space is: 



ipLtTfiidp + iA^ipL- / ipRd* Jdfi+iA^ipR (1) 



where ji, 



1, 2, (7i = 1 and a 2 = i. ipR,L,ipR,L 



are Grassmann variables and A^ is a U(l) gauge 
field. The model is defined on a torus of fixed 
physical size / x / with periodic boundary condi- 
tions for the electric field, E = e^d^A^, and the 
fermions. The topological charge q is an integer 
given by q = J E . The expectation value of 
the fermion bilinear (\?\?) is 



^[dipdipdA^e-^^^) 



N 
D 



(2) 



N has contributions only when q = ±1 and D has 
a contribution only when q = 0. 

The continuum torus is replaced by a toroidal 
square lattice with L sites in both directions. The 
lattice spacing a is l/L. Gauge fields U^ix) = 
e laA i*( x ) are associated with the links. In the over- 
lap formalism, N and D are expressed as follows: 



N 



-s g (u) 



,\L+)i 



D = J JJdU^x) 



e WML 



\L+)i 



(3) 



(4) 



* Speaker 



The a^ ^-term in N will differ from zero only for 
gauge fields carrying topological charge +1 and is 
the overlap expression for {ipRipL) while the a], a — 
term contributes only for charge —1 and repre- 
sents (iPl^r)- In our simulations we could have 
looked at only one of these terms; we took both 
in order to increase statistics. It is easy to prove 
that the two overlap matrix elements have exactly 
equal magnitudes. The CP violating angle 6 is set 
to zero. The states \L±)u are the ground states 
of many body Hamiltonians 

W ± = ^ a^ a H ± (a;,a,t/,/3;?7)aj /)/ 3 

x,a,y ,{3 

{4,a)%,/3j = 6a,t36 X ,y (5) 

a,/3 = 1,2, x = (xi,X2), = 0,1, ...,L— 1. 
The single particle hermitian hamiltonians W.^ 



2 



are given in [1]. We say that the topological 
charge q equals ±1 when the number of single 
particle states filled in and \L—)u differ 

by±l. 

The gauge action appearing in the formulae for 
N and D is a single plaquette action: 

S g (U) = 52s g (U p ), (6) 

p 

where U p is the oriented product of the four link 
elements making up the plaquette p. With such 
an action for a pure gauge theory the plaquettes 
become almost independent variables. To mini- 
mize the finite lattice spacing effects arising from 
S g (U), we chose the "heat-kernel" form [3]: 

e -s g (u p ) = ^ e-^[U p \ m (7) 

This choice assures exact scaling of the string ten- 
sion with the coupling g in the pure gauge theory. 
Other choices with similar properties are possible. 

Conventional MC techniques would generate a 
gauge field distribution corresponding to the in- 
tegrand in D (c.f. (4)) and thus compute the 
fermion bilinear in (2). But this will not work 
here because the set of gauge fields that con- 
tribute to the integral in (3) (viz., q = ±1) are 
completely disconnected from the set of gauge 
fields that contribute to the integral in (4) (viz., 
q = 0). Since the overlap formalism maintains 
this separation of gauge fields into topological 
classes even on a finite lattice we cannot use the 
conventional MC techniques. Of course, it is 
precisely the ability of the overlap formalism to 
maintain this property that enables us to measure 
the fermion bilinear on a finite lattice directly. 
We compute N/D in the following manner: Mea- 
sure N and D separately, as observables in a pure 
gauge theory, and calculate (\?\?) as 

<**) = § = f = §, Z= f[dU]e- s ^HS) 
z •* 

It is easy to generate statistically indepen- 
dent configurations for a single plaquette action: 
We choose a gauge and draw L 2 — 1 indepen- 
dent plaquette variables U p , each with probabil- 
ity e -^ (u "\ This fixes L 2 - 1 out of the L 2 + 1 



link variables. The last two link variables repre- 
sent Polyakov loops winding around the torus and 
since they do not affect the action we draw their 
values at random. The last plaquette is now fixed 
and comes with weight e - * 9 ^ 1 '"'*). The whole 
configuration is kept with probability e~ Sg ( Ul °- at \ 
Since the quantity we measure is strictly gauge in- 
variant on a finite lattice any choice of the gauge 
will give the same result. 

Because we treat the overlap factors as observ- 
ables in a pure gauge theory we must choose our 
parameters so that these factors do not fluctu- 
ate too rapidly. Since it has been shown [1] that 
for smooth gauge fields the overlap agrees with 
the continuum we use the continuum result of 
[4] to estimate the parameters of our simulation. 

We find that typical contributions are ~ e~ 'i^ 2 , 
where n = is the photon mass, and therefore, 
1 < ill < 3 is acceptable. 

The simulation is done in double precision. The 
single particle hermitian hamiltonians are di- 
agonalized using the Jacobi algorithm (scales as 
L 6 ). The determinants are computed using LU 
decomposition (also scales as L 6 ). The Jacobi 
procedure has a larger prefactor and computer 
time for the whole computation scales a little 
worse than L 6 for L k, 4 — 10 rather than the 
asymptotic L 8 (the last factor of L 2 comes be- 
cause of the x summation in (3)). 

We present our results in terms of the dimen- 
sionless quantity 

Results for nl = 3 have appeared in [2]. Within 
statistical errors (five percent) we found an agree- 
ment with the continuum result of [4] for L > 6. 

The error in [2] is dominated by the error in 
generating the correct distribution of the various 
topological sectors. We can isolate this contribu- 
tion to the error by first rewriting (8) as 

— N' N" 

= _ = R (10) 

x 1 D' D" v ' 

where N" is the same as N' but evaluated over 
the restricted set of gauge field configurations 
with q = ±1, and D" is the same as D 1 but 
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evaluated over the restricted set of gauge field 
configurations with q = 0. R is the ratio of the 
number of gauge field configurations with q = ±1 
divided by the number of gauge field configura- 
tions with q = 0. Until now we have been using 
the overlap definition of the topological charge. 
With this definition a separate evaluation of R 
would be costly. However, we have found that 
the geometrical definition q = J2 p — ir < <f> p < 
TT,U P = e l ^ v agrees with the overlap definition 
better than 99% for lattices with L > 6 if we re- 
duce the value of /il focusing on systems of smaller 
physical size, but with a larger ultraviolet cut- 
off. We can therefore evaluate R cheaply using 
the geometric definition or simply take it to have 
its continuum value from [4]. We choose the lat- 
ter approach and this reduces the errors signif- 
icantly. The ill = 1 results are the diamonds 
in fig. 1. In that figure the square represents 
the continuum result of [4]. Now the data devi- 
ates from the continuum prediction by about ten 
percent. The explanation for the deviation from 
the continuum result is simple and should have 
been anticipated: The "pure" Schwinger model 
does not include a marginal four- Fermi, Thirring, 
term. There is no symmetry that would turn this 
term off and any regularization will generate some 
unknown amount of it. Since it is marginal it 
does not vanish in the continuum and actually 
induces significant effects. The higher precision 
has allowed us to refine the conclusions of [2] and 
explicitly show the presence of a Thirring term in 
our regularization. 2 

The effect of the Thirring term is expected to 
be enhanced if one studies a model with more 
flavors than one. This can be done in the over- 
lap formalism at no extra cost: Simply raise each 
overlap factor that appears in eq. (3), (4) to a 
power given by the number of flavors Nf . We ob- 
tain results similar to the ones in fig. 1, and, in 
particular, the points do not approach the con- 
tinuum limit predicted by the exact solution of 
the continuum Schwinger model, providing fur- 
ther evidence for the presence of a Thirring term. 



The continuum Schwinger model with a 
Thirring term has been solved analytically in [5]. 
We leave the Thirring coupling constant, gx, 
and an overall normalization factor (reflecting the 
infinite renormalization when ^ 0) as free 
parameters and we perform a global fit of the 
continuum expression, modified for a finite lat- 
tice, to our data. We have collected data for 
L = 4,5,6,7,8,9,10, N s = 1,2,3,4 and /il = 
1, 1.5, 2, 3. The ill = 1, 1.5, L > 6 data for all N f 
are closer to the scaling region and were used for 
the global fit. The rest are used to estimate the 
"true" scaling violations (subleading in the lattice 
spacing). The solid line in fig. 1, shows the fit for 
N f = 1, ill = 1. The fit has a x 2 /d.o.f. = 1.09. 
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Figure 1. f L (l^l) 



2 As far as we know this issue never came up in previous 
lattice work concerning the Schwinger model. Maybe the 
accuracy was not good enough. 
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